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A SHORT NOTE ON A BERNSTEIN-BEZIER BASIS FOR THE 

PYRAMID 


JESSE CHAN, T. WARBURTON 

Abstract. We introduce a Bernstein-Bezier basis for the pyramid, whose restriction to the face 
reduces to the Bernstein-Bezier basis on the triangle or quadrilateral. The basis satisfies the standard 
positivity and partition of unity properties common to Bernstein polynomials, and spans the same 
space as non-polynomial pyramid bases in the literature ED EH m [ 9 ]. 

1. Introduction. Bernstein-Bezier bases have long been ubiquitous in graphics 
and computer-aided design muni, though they have recently received attention for 
their utility in the meshing of curved geometries and the numerical solution of partial 
differential equations mi nn mm]- Ainsworth and Kirby both noticed that the 
structure of Bernstein-Bezier polynomials on simplices allowed for linear-complextiy 
algorithms for the assembly and multiplication of finite element matrices Q3I2I2I3, 
as well as for linear-complexity Discontinuous Galerkin solvers [20j . The hierarchi¬ 
cal structure of Bernstein polynomials also facilitates the construction of structure¬ 
preserving vectorial basis functions HU 1. Additionally, their close connection to 
B-splines has been exploited to simplify the implementation of NURBS-based finite 
element methods [7j. 

Bernstein bases on tensor product elements (quadrilaterals and hexahedra) enjoy 
a natural tensor-product construction, while the construction of Bernstein bases on 
simplices rely on barycentric coordinates to generalize the construction in one space 
dimension. Bernstein-Bezier bases on the prism may similarly be expressed as the 
tensor product of triangular and one-dimensional basis functions. However, Bernstein- 
Bezier bases for the pyramid (which acts as a transitional piece to couple together 
hexahedral and tetrahedral elements 00) have received less attention. Recent work 
has considered the extension of such ideas to pyramids, but define Bernstein-Bezier 
pyramid bases (which are distinct from Bernstein-Bezier pyramid algorithms [Tit T]) 
by splitting the pyramid into two tetrahedra nnim. The work presented here presents 
an alternative construction based on pyramid bases found in the literature 0122 E3, 
which are non-polynomial, but contain the space of polynomials of total degree N 
under vertex-based mappings. Quadrature and the computation of matrices for finite 
element methods are also discussed, and condition numbers are presented for mass 
and stiffness matrices on the reference element. 

2. Bernstein-Bezier bases for triangles and quadrilaterals. Bernstein- 
Bezier basis functions on the triangle m may be defined in terms of the barycentric 
coordinates A 1; A 2 ,A 3 

tdN _ /~iN \i \j \k 

n ijk ~ ^ijk A 1 a 2 a 3 

with i + j + k = N and 

C N N\ 

ijk i\j\k\' 

For the unit right triangle with coordinates (r, s), these barycentric coordinates are 
given explicitly as 

Ai = 1 - (r + s), A 2 = r, A 3 = s. 
l 



Assuming the Duffy transform 


r = a(l — 6), s = b , 

which maps the unit cube with coordinates (a, b) £ [0,1] 2 to the unit right triangle, 
we have that 


Ai = (1 - a)(l — b), A 2 = a(l-6), A 3 = b. 


This gives a definition for the Bezier triangle basis on the quad 


B? jk = C$ k at(l-a) i bf'(l-b) i+ *. 


Since N = i + j + fc, i + j = N — fc, and we may rewrite the above as 
B? jk (a,b) = C» k ai(l-ayb k (l-bf- k 

= xy f ^ (1 - 

= Bf-yoXW = -Bf-*(oX(6) 

where 


B£(b)= u b\l-b) 


\N—k 


is the standard order N ID Bernstein polynomial and B^~ k {a) is the ith Bernstein 
polynomial of order N—k. This decomposition is at the heart of the sum-factorization 
techniques used in US El ED- Since the same barycentric approach is used to define 
Bernstein polynomials on a simplex of arbitrary dimension, the restriction of tetra¬ 
hedral Bernstein-Bezier basis functions to a triangular face results in the triangular 
Bernstein-Bezier basis. 

Bernstein-Bezier basis functions on quadrilaterals or hexahedra are defined using 
a tensor product construction; for example, on a reference quadrilateral with coordi¬ 
nates (a, b), 

B? j {a,b) = B?(a)B?(b). 

Similarly to the tetrahedron, restricting hexahedral Bernstein-Bezier functions to a 
single face results in Bernstein-Bezier polynomials over the quadrilateral. 

3. A Bernstein-Bezier basis for the pyramid. To extend this construction 
to the pyramid, we use a collapsed coordinate system m We define the unit cube 
with coordinates (a, 6, c) £ [0, l] 3 , such that the unit cube is mapped to the unit right 
pyramid through the transform 


r = a(l — c), s = b( 1 — c), t = c. 


A Bernstein-Bezier basis for the pyramid may be defined over the unit cube as follows 

B ijk (a,b,c) = B?-\a)Bf-\b)B% (c). 

where the indices obey 


0 < k < N, 0 < i, j < N — k. 
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The total dimension of this space is N p = (N + 1 )(N + 2)(2N + 3)/6. 

This construction is exactly what results from combining triangle Bernstein Bezier 
basis functions in the a, c and b, c coordinates. The Bernstein pyramid basis above 
satisfies the standard positivity and partition of unity properties, and the traces re¬ 
duce to Bernstein-Bezier basis functions on the triangle or quadrilateral. This latter 
property simplifies the enforcement of conformity between pyramids and tetrahedral 
or hexahedral elements. Additionally, the Bernstein pyramid space spans the same 
space as the rational basis of Bergot, Cohen, and Durufle 0, and thus contains poly¬ 
nomials of total degree N under a vertex-based mapping of the unit right pyramid to 
physical space. 

Lemma 3.1. The Bernstein pyramid basis satisfies a pointwise positivity and parti¬ 
tion of unity property on the pyramid. 

Proof. Since each component of the Bernstein pyramid basis is pointwise positive on 
the cube, the basis is pointwise positive over the pyramid. We may also show that 
the partition of unity property for the Bernstein pyramid is preserved: 

N N-kN-k N N-k N-k 

EEE = Ew E s E fc («) E = 1 

k =0 i—0 j—0 k =0 i—0 j =0 

by the fact that , B^ k , B™~ k all satisfy a partition of unity property over [0,1]. 

□ 

Lemma 3.2. The trace spaces of the Bernstein pyramid are Bernstein polynomials 
on the faces. 

Proof. We may prove the first property on the unit cube by restricting a, b to either 
0 or 1 for the triangular faces, and restricting c = 0 for the quadrilateral face. For 
the quadrilateral face, since Bk(c) = 0 for k > 0, the Bernstein basis is nonzero only 
if k = 0. 


Bij o = B?{a)B*f{b) 

which is exactly the tensor product Bernstein basis over the quadrilateral face. For 
the triangular faces, we take a,b = 0,1. We show the trace space for a = 0; the other 
faces are similar. If a = 0, the Bernstein basis is nonzero only for i = 0, and the 
nonzero Bernstein basis functions follow the form 

B 0j k = Bf~ k {b)B»{c). 

This is identical to the Bernstein basis on the triangle after a mapping to the unit 
quadrilateral. □ 

Lemma 3.3. The Bernstein pyramid space is identical to the pyramid space of Bergot, 
Cohen, Durufle W- 

Proof. This is simplest to show by showing equivalence with the pyramid basis pre¬ 
sented in [5], which is also equivalent to the basis of Bergot, Cohen, and Durufle. 

Ma, b , c) = £ k (a)£ k (b) ^(c), 
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or equivalently 


/, \N-k 

e?- k (a)e?- k (b) pr~ k)+3 (c). 

As the spaces spanned by each basis are of the same dimension, it simply remains to 
show that their span is identical. Since the functions in a, b 

B?- k {a)Bf~ k {b), e?- k (a)e?~ k (b), 0 <i,j<N-k 

both span Q^-kia, b), and the functions in c 

(—) P 2 k {N - k)+3 (c), B k (c) 

are both linearly independent homogeneous polynomials of total order N, the two 
bases span the same approximation space. 

□ 


As a result, the Bernstein pyramid contains the space of polynomials P N on any 
vertex-mapped pyramid, and high order accurate approximations may be constructed 
on general vertex-mapped pyramids. 

Unlike the orthogonal bases for pyramids constructed previously, the Bernstein 
basis absorbs the extra homogenizing factors of c k or (1 — c) k into the definition of 
the Bernstein-Bezier basis in the c direction, resulting in a very concise formula. 

3.1. Derivatives, quadrature, and evaluation of matrices. To compute 
derivatives with respect to reference coordinates r, s, t , we use the chain rule, involving 
factors of 


da db 1 da r a 

dr ds 1 — c dt (1 — t) 2 1 — c 

db _ s b 

dt (1 — t) 2 1 — c 

If k = N, the derivative is zero. For k < N, we have 


dB?-*(a) 

da 


B 


dB%k 

dr 


9B m _ B N- k(n] dB > 
~dT- Bi {a) — 


(b) 


BUc) 


N-k 

j v ~' 1 — C 
N-k/ 


\b) B%{c) 


db 


1 — c 


9B? jk dB% k 

( a \ 

I + dBijk 1 

( b ) 


dt da ' 

U -c) 

l+ db 1 

U-cJ 

l+ dc 


The additional factors of a,b and 1/(1 — c) may also be absorbed into the Bernstein 
polynomials, and properties of one-dimensional Bernstein-Bezier polynomials may be 
used to rewrite all derivatives in terms of Bernstein polynomials. However, doing so 
does not readily simplify expressions for derivatives, which do not appear to be as 
sparse or concise for the pyramid as they are for the tetrahedron. 

Under the collapsed coordinate system, integrals on a physical pyramid V are 
computed via the transformation 


u dx d y d 2 



c) 2 da d& dc. 
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where J is the determinant of the Jacobian of the mapping from reference pyramid 
V to V. For vertex-mapped pyramids, J, as well as all change of variables factors 
are bilinear in the a,b coordinates and constant in c [SJ Lemma 3.5]. An 
appropriate quadrature for the pyramid may then be constructed using the tensor 
product of one-dimensional Gaussian quadratures in the a, b coordinates and Gauss- 
Jacobi quadratures with weights (2,0) in the c coordinate, which integrates exactly 
the above expression. 

The entries of the mass matrix 


Mijk,Lmn — / ^ijk^lmn — / 4*ijk ( t ) lmnJ 
J-p Jv 

= [ B?- k {a)B»- k {a) f Bf~ k {b)B^~ k {b) f B»(c)B%(c)(l - cfj 

J a J b J c 

are straightforward to compute. Fast algorithms for Bernstein-Bezier bases [I8J GHEE] 
which use that the product of two Bernstein polynomials is again a scaled Bernstein 
polynomial, and that the integrals of Bernstein polynomials are explicitly known 
in terms of binomial coefficients — may also be applied to the evaluation of the 
pyramidal Bernstein-Bezier mass matrix, though unlike vertex-mapped tetrahedron, 
the non-constant factor J must be taken into account. The weak derivative matrices 

s x ,sy 7 s z 


QX 


oV 

^ijk,lmn 


C2 

u ijk,lmn 


l 

l 

i 


tdN un lmn 

ijk dx 

tdN UJD lmn 

ijk dy 

r) 

r?N un lmn 

ijk dz 


are also computed exactly using the quadrature rule described above for vertex- 
mapped pyramids 0. However, the stiffness matrix is not integrated exactly using 
quadrature unless the transformation is affine (i.e. the base is a planar parallelogram) 
due to rational factors in the resulting integrals 0. A non-conforming error analysis 
and numerical experiments appear to indicate that high order accuracy may still be 
attained using the quadrature rule described above to compute the stiffness matrix 
0- 


3.2. Conditioning of matrices. Figure 3.1 reports computed condition num¬ 
bers of the reference mass and stiffness matrices for both pyramidal and tetrahedral 
Bernstein-Bezier bases at various orders of approximation. Dirichlet boundary con¬ 
ditions are enforced to ensure that the stiffness matrix is nonsingular. The condition 
number of each matrix grows exponentially in the order N, as expected for Bernstein- 
Bezier type bases. However, the growth of the condition numbers for the pyramid 
is more rapid than the growth of the condition number for the tetrahedra, and may 
need to be addressed in order to maintain numerical accuracy at high orders of ap¬ 
proximation. 
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Fig. 3.1: Condition numbers of reference tetrahedral and pyramidal Bernstein-Bezier 
mass and stiffness matrices at various orders of approximation N. 
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